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Abstract 

We study a recurrence defined on a three dimensional lattice and prove that its 
values are Laurent polynomials in the initial conditions with all coefficients equal to 
one. This recurrence was studied by Propp and by Fomin and Zelivinsky. Fomin and 
Zelivinsky were able to prove Laurentness and conjectured that the coefficients were 
1. Our proof establishes a bijection between the terms of the Laurent polynomial and 
the perfect matchings of certain graphs, generalizing the theory of Aztec diamonds. In 
particular, this shows that the coefficients of this polynomial, and polynomials obtained 
by specializing its variables, are positive, a conjecture of Fomin and Zelevinsky. 

1 Historical Introduction 

The octahedron recurrence is the product of three chains of research seeking a common 
generalization. The first is the study of the algebraic relations between the connected minors 
of a matrix, and particularly of a recurrence relating them known and Dodgson condensation. 
Attempting to understand the combinatorics of Dodgson condensation lead to the discovery 
of alternating sign matrices and Aztec diamonds. Aztec diamonds are graphs whose perfect 
matchings have extremely structured combinatorics and soon formed their own, second line 
of research as other graphs were discovered with the similar regularities. The third is the 
study of Somos-sequences and the Laurent phenomenon, which began with an attempt to 
understand theta functions from a combinatorial perspective. In this section we will sketch 
both lines of research. In the following section, we will begin to describe the vocabulary and 
main results of this paper. 
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1.1 Dodgson Condensation and the Octahedron Recurrence 

Let rriij be a collection of formal variables indexed by ordered pairs G Z^. Let ii, i2, ji 
and j2 be integers with i2 — h = 32 — ji > and set 

D^f = det(mi,). .... . .... . 



Charles Dodgson Dodg| observed that 



7-)i2j2 7-)(i2 + l),02 + l) _ 7-)(«2 + l),02 + l) p,i2,j2 _ 7-)«2,(i2 + l) r)(22+l)j2 

suggested using this recursion to compute determinants and found generalizations that could 
be used when the above recursion requires dividing by 0. (This formula is not original to 
Dodgson, according to page 7 of |Lunj this formula has been variously attributed to Sylvester, 



Jacobi, Desnanot, Dodgson and Frobenius. However, the use of this identity iteratively to 
compute connected minors and the detailed study of this method appears to have begun 
with Dodgson and is known as Dodgson condensation.) 

Before passing on, we should note that the study of algebraic relations between determi- 
nants and how they are effected by various vanishing conditions is essentially the study of 
the flag manifold and Schubert varieties. This study lead to the invention of cluster algebras 
which later returned in the Laurentness proofs we will discuss in Section 11.31 See |FomZel2] 
for the definitions of Cluster Algebras. The precise relation between cluster algebras and 
(double) Schubert varieties is discussed in |BFZj . although the reader may wish to first 
consult the references in that paper for the history of results that preceeded it. 

One can visualize this recurrence as taking place on a three dimensional lattice. The 
values rriin are written on a horizontal two dimensional plane. The values D^^i^ live on a 
larger three dimensional lattice with -D*^j^ written above the center of the matrix of which 
Z)-JjJ is the determinant and at a height proportional to 12 — ii = — ji- The six entries 
involved in computing the recurrence lie at the vertices of an octahedron. We make a change 
of the indexing variables of the lattice so that the lattice consists of those {n,i,j) with 
n + i + j = mod 2. Then the recurrence above can be rewritten as 

f{n,i,j)f{n-2,i,j) = f{n-l,i-l,j)f{n-l,i + l,j)-f{n-l,i,j - l)/(n - 1, i, j + 1) 

with /(O, = rrii+itd. and /(—I, = 1. 

2 2 

Robbins and Rumsey [RobRum] considered modifying the recurrence to 

f{n,i,j)f{n-2,i,j) = /(n - 1, i - 1, j)/(n - 1, i + 1, j) + A/(n - 1, i, j - l)/(n - 1, i, j + 1) 

and also generalized by allowing /(— j) to be arbitrary instead of being forced to be 1. 
We write x{i,j) = f{0,i,j) when i + j is even and x{i,j) = f{—l,i,j) when i + j is odd. 
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Robbins and Rumsey discovered and proved that, despite these modifications, the functions 
f{n,i,j) were still Laurent polynomials in the x{i,j) and in A, with all coefficients equal to 
1. 

Robbins and Rumsey discovered that the exponents of the x{i,j) in any monomial occur- 
ing in a /(no, io,jo) formed pairs of "compatible alternating sign matrices." We will omit the 
definition of these obects in favor of discussing a simpler combinatorial description, found 
by Elkies, Kuperberg, Larsen and Propp |EKLPj . in terms of tilings of Aztec diamonds. 



1.2 Aztec Diamonds 

If G is a graph then a matching of G is defined to be a collection M of edges of G such 
that each vertex of G lies on exactly one edge of M. Matchings are also sometimes known 
as "perfect matchings", "dimer covers" and "1-factors". Fix integers no, io and jo with 
n-o + ^0 + Jo — mod 2. In this section, we will describe a combinatorial interpretation of 
the entries of f{no,io, jo) in terms of matchings of certain graphs called Aztec diamonds. 
For simplicity, we take A = 1 in the recurrence for /. 

Consider the infinite square grid graph: its vertices are at the points of the form [i + 
1/2, j + 1/2) with (z,j) e Z^, its edges join those vertices that differ by (0, ±1) or (±1,0) and 
its faces are unit squares centered at each G Z^. We will refer to the face centered at 
(i,j) as («, j). Consider the set 5'(no,io, jo) of all faces (i, j) with \i — io\ + \ j — jo\ < uq. Let 
G(no, io,jo) be the induced subgraph of the square grid graph whose vertices are adjacent to 
some face of >S'(no, io, jo)- We will call this the Aztec diamond of order no centered at (io, jo)- 

We define the faces of G(no,'io, jo) to be the squares (i, j) for |i — iol + |j ~ jo| ^ n-o- Note 
that this includes more squares than those in S'(no, io, jo); in particular, it contains some 
squares only part of whose boundary lies in G{no, io, jo)- We refer to the faces in ^(no, io, jo) 
as closed faces and the others as open faces. Let M be a matching of G'(n,o, io, jo)- If (i, j) 
is a closed face then either 0, 1 or 2 edges of (i, j) are used in M; we define e(i, j) to be — 1, 
or 1 respectively. If x{i,j) is an open face then either or 1 edges will be used and, we 
define e(i, j) to be or 1 respectively. Set m(M) = Yl^ih jY^^'''\ where the product is over 
all faces (i, j) of (^(no, io, jo)- 

The Aztec Diamond Theorem, /(no, io, jo) = X] where the sum is over all match- 

ings M o/G(no,io, jo)- 

Proof. It seems to be difficult to find a reference which states this result exactly in this form. 
IRobRurn] shows that /(no,io,jo) is equal to a sum over "compatible pairs of alternating 
sign matrices" and |EKLPj shows that such pairs are in bijection with perfect matchings 



of the Aztec Diamond. Tracing through these bijections easily gives the claimed theorem. 
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Another proof can be found by mimicing the method used in |Kuoj to count the matchings 
of the Aztec Diamond. This result will also be a corollary of this paper's main result. □ 

Corollary 1. The order n Aztec diamond has 2"("+^)/^ perfect matchings. 

Proof. Clearly, if we plug in 1 for each x{i,j), then m{M) = 1 for every matching M so the 
number of matchings of the order n Aztec diamond will be f{n,i,j)\x(^ij^=i. We will denote 
this number by g{n). Then the octahedron recurrence gives g{n)g{n — 2) = 2g{n — 1)^ and 
the claim follows by induction. □ 

Many other families of graphs have been discovered for which the number of matchings of 
the n*^ graph is proportional to a constant raised to a quadratic in n. They are summarized 
and unified in |Ciucuj . One of the purposes of this paper is to give proofs of these formulae 
as simple as the corrolary above. 

The monomial m(M) keeps track of the matching M in a rather cryptic manner; it 
is not even clear that m(M) determines M, although this will follow from a later result 
(Proposition [T7j). It is natural to add new variables that code directly for the presence or 
absence of individual edges. Let i and j be integers such that i + j = 1 mod 2. Draw the 
plane such that i increases in the east-ward direction and j in the north-ward. Let the edges 
east, north, west and south of (i, j) be labelled a{i,j), b{i,j), c{i,j) and d{i,j) respectively. 
Introduce new free variables with the same names as the edges. 

We now define an enhanced version of / by 

f{no,io,jo) 

M eGM 

where the sum is still over all matchings of the order no Aztec diamond. Each matching is 
now counted by both the monomial m{M) and the product of the edges occuring in that 
matching. The new / obeys the recurrence 

fin,i,j) = {a{i + n-l,j)c{i-n + l,j)f{n-l,i,j + l)f{n-l,i,j -1) + 

b{t,j +n- l)dit,j -n+ l)f{n - 1, z + l,j)f{n - 1, z - 1, j)) /f{n - 2, z, j) 

1.3 Somos Sequences and the Octahedron Recurrence 

Somos has attempted to rebuild the theory of Theta functions on a combinatorial basis. 
Theta functions are functions of two variables, traditionally called z and q. The simplest 
Theta function is Q{z,q) = Xl^-oo • Fro^i now on, we will concentrate only on 

the dependence on z and so write Q{z). For a general introduction to 6 functions, see for 
example chapter XXI of |Whit Wat] . 



4 



Somos began with a well known result, that, for any xq and Zi, Q{z) would obey 



Q{zo + nzi)Q{zo + {n - 4:)zi) = rQ{zo + {n - 8)2:1)6(^0 + {n - l)zi) + sQ{zo + {n - 2)zi 



i2 



for r and s certain constants depending in a complicated manner on q. Somos proposed using 
this recurrence to compute Q{zo + nzi) in terms of r, s and Q{zo), Q{zo + zi), 9(2;o + 22;i) and 
6(2:0 + 82:1). He discovered that 0(2:0 + nzi) was a Laurent polynomial in these 6 variables. 
This is not difficult to prove by elementary means, but suggests a deeper explanation may be 
lurking. Moreover, he discovered that all of the coefficients of this Laurent polynomial were 
positive but was unable to prove this; this will be proven for the first time as a corollary of 
the results of this paper. 

After Somos's work, many other recurrence with surprising Laurentness properties were 
discovered; see |Galej for a summary. The family of most importance for this paper is the 
Three Term Gale- Robinson Recurrence: Let k, a and b be positive integers with a,b < k. 
The Three Term Gale-Robinson Recurrence (abbreviated to "Gale-Robinson Recurrence" in 
the remainder of this paper) is 

g{n)g{n — k) = rg{n — a)g{n — k + a) + sg{n — b)g{n — k + b). 

Note that the Somos recurrence is the case {k, a, b) = (4, 1, 2). Note that if {k, a, b) = (2, 1, 1) 
and r = s = 1, we get the recurrence g{n)g{n — 2) = 2g{n — l)'^, the recurrence for the number 
of perfect matchings of the Aztec diamond. 

The Gale- Robinson recurrence, for every value of {k, a, b), can be thought of as a special 
case of the Octahedron recurrence. This observation was first made by Propp and seems to 
have first appeared in print in [FomZell] . The reduction is performed as follows: suppose that 



g obeys the Gale-Robinson recurrence. Define f{n,i,j) by f{n,i,j) = (^(^ ^"+(^° fc>+(2b k)j -.^ 
Then / obeys 

fin,ij)fin-2,i,j) = rf{n-l,i-l,j)f{n-l,i + l,j) + sf{n-l,i,j-l)f{n-l,i,j + l). 

Thus, if we choose a{i,j), b{i,j), c{i,j) and d{i,j) to be constants a, b, c and d with 
ac = r and bd = s, f obeys the octahedron recurrence. In these circumstances, we must 
study the octahedron recurrence not with the initial conditions f{0,i,j) and /(— l,i,j) but 
with the initial coniditions f{n,i,j) for {n,i,j) such that —k < < 

It therefore beomes our goal to study the octahedron recurrence with general initial 
conditions. The main result of this paper will be that, with any inital conditions, f{n, i,j) is 
a Laurent polynomial in the initial coniditions and in the coefficients a{i,j), b{i,j), c{i,j) and 
d{i,j). Moreover, all of the coefficients of this Laurent Polynomial are 1. In particular, they 
are positive. In the reduction of the Gale-Robinson recurrence to the octahedron recurrence. 
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we imposed that many of the initial conditions of the octahedron recurrence be equal to each 
other. Thus, many monomials that would be distinct if the octahedron recurrence were run 
with arbitrary inital conditions become the same in the Gale-Robinson recurrence. However, 
we can still deduce that all of the coefficients of the Laurent polynomials computed by the 
Gale-Robinson recurrence are positive. 

2 Introduction and Terminology 

In this section we will introduce the main objects and results of this paper. The main object 
of study of this paper is a function / defined on a three dimensional lattice by a certain 
recurrence. We denote the set of initial conditions from which / is generated by X and 
consider varying the shape of X inside this lattice. It turns out that the values of / are 
Laurent polynomials in the initial variables where every term has coefficient 1. We are able 
further more to find families of graphs such that / gives generating functions for the perfect 
matchings of these graphs. Special cases of this result include being able to choose X so 
that / gives us the values of three term Gale-Robinson sequences or such that the families 
of graphs are Aztec Diamonds, Fortresses and other common families of graphs with simple 
formulas for the number of their matchings. 

In the next subsection, we introduce the vocabulary necessary to define /. In the following 
subsection, we introduce the vocabulary necessary to define the families of graphs. 

2.1 The Recurrence 

Set 

£ = {(n, i, j) e Z^, n — i^2 mod 2} 

and set 

^ = Q') e X {«,^c, (i}, i + j = 1 mod 2} 

T = {{t,j)eZ^} 

where q denotes one of the four symbols a, b, c and d. C, S and T stand for "lattice", 
"edges" and "faces." 

We call and (^2,^2) £ ^ "lattice- adjacent" if |(ii — ^2)! + |(ji —^2)1 = 1- (Another 

notion of adjacency will arise later.) 
For {no,io,jo) e >C, set 

P{no,io,jo){hj) =no-\i-io\-\j - Jo\. 
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(We will often drop the subscript on p when it is clear from context.) Let 

C(„o,io,jo) = {{n,i,j) e C : n < P(no,io,jo)iiJ)} 

C{no,io,jo) = {{n,i,j) e C : n < P(no,io,jo)ihj)} 

dC = {{n,i,j) e C:n = p^no,io,jQ)ihj)} 

We call C, C and dC the cone, inner cone and outer cone of (no,«o, jo) respectively. C is a 
square pyramid with it's vertex at (no,io,io)- 
Let h: J-" Z and define 

1 = {(^,J» e C, n = h{i,j)} 
W = e £, n > h{i,j)}. 

Assume that 

1. \h{ii,ji) — h{i2,j2) \ = 1 if (^1, ji) and {12,32) are lattice-adjacent. 

2. {h{i,3),i,3)eC. 

3. lim|i|+|j|^oo/i(^, j) + \A + lil = oo- 

The last condition is equivalent to "For any {n,i,j) G C, C(^n,i,j) H W is finite." We will 
call such an h a "height function" and call a function h that obeys the first two conditions 
a pseudo-height function. Note: this definition is not related to the height functions in the 
theory of Aztec diamonds, described in |CEPj . 

Let K be the field of formal rational functions in the following infinite families of variables: 
the family x{i,j), i, 3 G I? and the families a{i,j), b{i,j), c{i,j), d{i,j) where i, 3 G I? 
and i + j = mod 2. (Clearly, JF indexes the x's and £ indexes the a's, 6's, c's and rf's). 
Let R <Z K he the ring Z[x(i, j), l/x(i, j), a{i,j), b{i,j), c{i,j), d{i,j)]. For reasons to appear 
later, we call the x's the "face variables" and the a's, 6's, c's and li's the "edge variables." 

We define a function f : 1 UU ^ K recursively by f{n,i,j) = x{i,j) for {n,i,j) G X and 

fin,i,j) = {a{i + n-l,j)c{i-n + lj)f{n-l,i,j + l)f{n-l,i,j -1) + 

b{i,3 +n- l)d{i,3 -n + l)/(n - 1,^ + 1,3) f{n - 1,2 - 1, j)) /f{n - 2,z,j) 

for {n,i,j) G U. Our third condition on h ensures that this recurrence will terminate. 
The main result of this paper is 

The Main Theorem (First Statement). f{n,i,j) G R and, when written as a Laurent 
polynomial, each term appears with coefficient 1. Moreover, the exponent of each face variable 
is between —1 and 3 (inclusive) and the exponent of each edge variable zs or 1. 

To describe our result more precisely, we need to introduce some graph-theoretic termi- 
nology. 
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2.2 Graphs With Open Faces 

Let Go be a connected planar graph with a specified embedding in the plane and let ei, 62, 
. . . e„ be the edges of the outer face in cychc order. We define a "graph with open faces" to 
be such a graph Go with a given partition of the cycle ei, . . . e„ into edge disjoint paths e^, 

^i+l ; ■ ■ ■ Gj . 

Denote by G the graph with open faces associated to Go and the partition above. We 
call a path e^, ■ ■ - Cj in the given partition an "open face" of G. We denote the open faces 
of G by -Fo(G). We refer an interior face of Go as a "closed face" of G and denote the set 
of them by Fc(G). We set F(G) = Fo{G) U Fc{G) and call a member of F(G) a face of G. 
Note that the exterior face of Go is not usually a face of G. 

The image associated to a graph with open faces is that open disks have had portions 
of their boundary glued to the outside of Go along the paths e^, ■ ■ - Cj with some additional 
boundary left hanging off into space. 

We will at times need to allow Go to be infinite, in which case Go might have several outer 
faces, or none at all. We will not describe the appropriate modifications, as they should be 
obvious. The infinite graphs which will come up in this paper will have no outer face and 
the corresponding graphs with open faces will have no open faces. 

We will refer to an edge or vertex of Go as an edge or vertex (respectively) of G. We say 
that an edge e borders an open face / of G if e lies in the associated path e,;, . . .ej. All other 
incidence terminology (cndpoint, adjacency of faces, etc.) should be intuitive by analogy to 
ordinary planar graphs. We use E{G) and V{G) to denote the edges and vertices of G. 

By a map G' ^ G we mean a triple of injections V{G') ^ V^(G), -E(G') "-^ E{G) and 
F(G') ^ F{G) compatible with the adjacency relations. We say G' is a sub-graph with 
open faces if such a map exists. Note that we do not require that open faces be taken to 
open faces. 

Let -R(G) be the ring Z[xf, l/xf,ye] where xj are formal variables indexed by / e F{G) 
and Ue are formal variables indexed by e e E(G). 

For the rest of this subsection, assume that G is finite. A matching of G is a collection 
of edges M such that every vertex lies on exactly one edge of M. 

Let M be a matching of G. For e G E{G), set S{e) = 1 if e G M, otherwise. For 
/ e F{G), let a be the number of edges of / that lie M and b the number of edges of / not 
in M. If / G F,(G), set 

if / G Fo{G) set 

, . \b-a' 
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(Note that e(/) > -1.) 

For any matching M, set 

m{M) =Y[y^^^''^Y[xf\ 

e f 

For any graph with open faces G, set 

m(G) =^m(M) 

M 

where the sum is over all matchings of G. We call m{M) the matching monomial of M and 
m{G) the matching polynomial of G. 

It will turn out that the actual graphs to which we will apply this definition are bipartite, 
so every closed face has an even number of edges and we may omit the [ ] in this case. 
However, it will be convenient to have a definition that is valid for all graphs. 

We can now give our second statement of the Main Theorem. 

The Main Theorem (Second Statement). For any height function h we can find an 
infinite graph with open faces Q , a decomposition E{Q) — Ey, U Eu and injective maps 
a: Eyj £ , a: F ^ T and a family of sub- graphs with open faces Gn,i,j indexed by {n, i,j) e 
U such that 

1. If we define a: R{Q) R by Q:{xf) = x{a{f)) for f G F{Q), a(|/e) = where 
a{e) — {i,j, q) for e e Ey^ and a{ye) — 1 for e e E^, then 

f{n,i,j) = a{m{Gn,i,j)). 

(We have implicitly used the obvious injection R{Gn,ij) ^ R{G)-) 

2. If {n',i',j') G Cn^ij then Gn'^i'j' is a sub-graph with open faces of Gn^ij 

We refer to the edges in Eyj as the "weighted edges" and the edges in E^ as the "unweighted 
edges". 

Clearly, this will imply the previous statement of the main theorem, except for the bounds 

on the exponents and the claim that the coefficients are 1. The first will follow by showing 
as well that every face of G has < 8 edges. The second will follow by showing that the edges 
of Eu are vertex disjoint, so a matching M is uniquely determined by M n E^j. 
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2.3 Plan of the Paper 

In Section El we will describe an algorithm we refer to as "the method of crosses and 
wrenches" for finding the graph Q and the subgraphs G(^n,i,j)- We will postpone proving 
the correctness of the algorithm to present applications of the theory in Section 13.51 In 
particular, we will carry out the examples of Somos-4 and Somos-5 and show explicitly the 
families of graphs associated to them. We will also show that, by choosing certain periodic 
functions for h, we can create fortress graphs and families of graphs studied by Chris Dou- 
glass and Matt Blum. The number of matchings of these graphs are powers of 5, 2 and 
3 respectively, in each case we will give a rapid proof of this by induction from our Main 
Theorem. Fortresses and Douglass' graphs are discussed in jCiucuj and these formulas are 
also proved there. 

In the Sectional we will first relate the matchings of the graphs G(^n.i.j) to the matchings 
of certain infinite graphs subject to "boundary conditions at infinity". This will remove 
the elegant property that G(n',i',j') ^ ^(n,i,j) when {n',i',j') G Ci^n,i,j), but it will create 
objects better suited to an inductive proof. We will then prove the main theorem by vary- 
ing I and holding {n,i,j) fixed. In the process, we will recover variants of the "Urban 
Renewal" operations of [Ciucui . 

In Section we will give a second proof, inspired by a proof of |Kuoj for the Aztec 
Diamond case, that works by holding X fixed and varying {n,i,j). This proof has some 
additional consequences that the first proof does not. However, due to the extraordinary 
number of cases that would otherwise be involved, we will not check that all of the exponents 
work out. 

In Section El we make some final comments. 

2.4 A Note on Simplification of Results 

One might wonder whether the introduction of so many families of variables is necessary to 
the paper. For most parts of the paper, one may think of any family variables one does not 
want to deal with as simply being set equal to 1 without difficulty. However, there is one 
place where this does not work: our first proof of the main theorem is inductive and the 
induction only goes through if the theorem is stated with the entire family of face variables 
present. 
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3 The Method of "Crosses and Wrenches" 

In this section we describe how to find the graphs Q and G(n,i,j) discussed above. As a 
running example, we we'll compute /(3, 1, 0) with h{i,j) = The goal is to predict the 

formula: 

/(3,1,0) = a(3,0)c(-l,0)a(2,-l)c(0,-l)x(l,-2)x(l,-l)~^x(l,l) 

+ a(3, 0)c(-l, 0)6(1, 0)rf(l, -2)x(l, 0)"^x(0, -l)a;(2, -l)x(l, -l)-^a;(l, 1) 
+ 6(1, 2)d(l, -2)a(l, 0)c(-l, 0)x(0, -l)a;(0, l)x(0, 0)"^x(2, 0)x(l, 0)"^ 
+ 6(1, 2)rf(l, -2)6(0, l)d{0, -l)x(-l, 0)x(0, 0)"^x(2, 0). 

The four terms in the above formula will eventually be shown to correspond to the four 
matchings of the graph with open faces shown in Figure ^ Here the dashed edges seperate 
the various open faces and the graph is drawn twice so that the face and edge labels will not 
overlap each other. The four matchings are shown in Figure |21 the numbers in Figure |21 are 
the exponents of the corresponding face variables. 

The octahedron recurrence with these initial conditions appears section 9.4 of ^FGj . A 
"tropicalized" version of our running example appears in jKTWj . 

3.1 The Infinite Graph Q 

Let h{i,j) be a height function. We describe the graph Q by describing its dual. The faces 
of Q (which are all closed) are indexed by the elements of I and the map a sends the face 
{n,i,j) G X to (z,j) G JF. Picture the face {n,i,j) as centered at the point (z,j) G M?. 

If (ii,ji) and («2?j2) are lattice-adjacent (i.e. —12) — (ji — i2)\ = 1) then (ni,zi,j2) 
borders {n2,i2,j2)- In addition, if \ii — i2\ = \ji — J2I = 1 then (ni,zi,ji) borders (n2,^2?j2) 
iff h{ii,j2) 7^ h{i2,ji)- No other pairs of faces border and faces which border border only 
along a single edge. 

We refer to this method of finding Q as the "method of crosses and wrenches" because it 
can be describe geometrically by the following procedure: any quadruple of values 

( h{t,j) h{z + l,j) \ 
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Figure 1: The Graph of our Running Example 
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must be of one of the following six types: 



f h h + l\ fh + 1 
\h + l h J \ h 

fh + 2 h + l\ fh + 1 
\h + l h J \h + 2 

In the center of these squares, we draw a 
fourth, and fifth and sixth cases respectively. 



h \ f h h + l\ 
h + l) yh + l h + 2) 

h \ fh + l h + 2\ 
h + l) \ h h+l) 

^, ^ror +^ in the first and second, third and 
'he diagram below displays the possible cases. 



h h + l\ (h+l h \ / h h + 1^ 

^h + 1 h j \ h h + l) \h + l h + 2, 

'h + 2 h+l\ /h + l h \ /h + l h + 2'' 

) ( ^ ) ( ^ 

yh + 1 h / \h + 2 h + l/ \ h h + l, 

We then connect the four points protruding from these symbols by horizontal and vertical 
edges. (We often have to bend the edges slightly to make this work. Kinks introduced in 

this way are not meant to be vertices, all the vertices come from the center of a ^ or from 
the two vertices at the center of a ^ or ^ .) We refer to the ^'s as "crosses" and the ^'s 
and ^'s as "wrenches" . 

In our running example, if i + j = then /^(^/j+i) h{i+i\'+i) ~ 10 place a 

Otherwise, we get a the resulting infinite graph, shown in Figure El consists of a 

diagonal row of quadrilaterals seperating a plane of hexagons. 

3.2 Labeling the Edges 

We will now describe the decomposition E = U Eu and the map a: E^ — > £. 

The set Ey^ will consist of the horizontal and vertical edges, i.e., those separating lattice- 
adjacent faces. The set E^ will consist of the diagonal edges, i.e., those which come from 
the center of a wrench. 

Consider any edge e of i?^. Such an edge lies between two faces (rii, h, ji) and (^2, ^2? j2) 
with \ii —i2\ + \ji — 72! = 1- Without loss of generality, {ii — 12) + (ji —32) = 1- (Otherwise, 
switch the names of (ii, ji) and {12^32)-) There are four cases: 
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1. If ii > 12 and rii > n2 then a{e) = {ii + n2,ji, a) = {12 + 'ni,j2, a). 

2. If ii > i2 and rii < n2 then a(e) = {ii —n2,ji,c) = {12 —ni,j2,c). 

3. If ji > j2 and ni > n2 then a{e) = + n2, b) = {12-, 32 + ni, h). 

4. If ii > j2 and rii < n2 then a{e) = — n2, a) = {12^32 — ni, d). 
The reader may wish to refer to Figure ^ 

3.3 Finding the Subgraphs G(^n^ij-^ 

Finally, we must describe the sub-graph with open faces G(n,i,j) of Q that corresponds to a 
particular {n,i,j). We will abbreviate G(^n,i,j) by G in this paragraph. The closed faces of G 
will be C(n,i,j) n X. The edges of G will be the edges of Q adjacent to some face in X fl G. 
The open faces of G will be those members of X some but not all of whose edges are edges of 
G. Note that every open face of G lies in dC(^n,i,j) H ^ but the converse does not necessarily 
hold. Note also that there are never any edges that separate two open faces; even if those 
two open faces are lattice-adjacent. 

Another way to describe the faces of G which is more intuitive but harder to write down 
is that the faces of G correspond to those f{n',i',j') which are used in compuing f{n,i,j) 
and the closed faces correspond to those f{n',i',j') which are divided by in the course of 
this computation. 

In our running example, the closed faces are (0,0,0), (1, 1,0) and (0, 1, —1). The open 
faces are (1, 0, 1), (2, 1, 1), (1, -1, 0), (2, 2, 0), (1, 0, -1), (1, 2, -1) and x(l, 1, -2). 

3.4 The Main Theorem 

We can now give the final statement of our Main Theorem. 

The Main Theorem. For any height function h, define Q , a and Gi^n.i.j) by the algorithm 
of the previous sections. Then the second statement of the main theorem holds with regard 
to these choices. 

3.5 Some Basic Facts 

It is easy to check that every closed face of Q has 4, 6 or 8 sides. (Simply check all possible 
values for the face's eight neighbors.) Hence all of the G(^n,i,j) are bipartite. Moreover, every 
face has < 8 edges, as promised above. 
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It is also easy to fulfill another promise and check that the unweighted edges are vertex 
disjoint: they lie in the middle of wrenches and do not border each other. 

It is clear that a: F{Q) — > is injective (and in fact, bijective). We now show: 

Proposition 2. The map a: £ is injective. 

Proof. Consider an edge e = (io, Jo? (^) ^ the cases of 6, c and d are extremely similar. Any 
edge with label e must be between two faces of the form (n, iQ—n — l, jo) and (n+1, iQ—n^jo). 
We must show there is at most one value of n for which (ra, io — — 1, Jo) and (n + 1, — Jo) 
both lie in 1. 

Suppose, for contradiction, there are two such values: n and n' . Without loss of generality, 
suppose that n < n' . Then 

h{iQ - n', jo) - K^o - n - I, jo) = {n' + I) - n > n' - n - I = {iQ - n') ~ {iq - n - 1). 

But h{i,j) must change by ±1 when i or j changes by 1, so h can not change by this 
much between (io — n',jo) and {io — n — 1, jo), a contradiction. □ 

We now describe where the faces of G(no,iojo) he in the plane. We make use of the 
following abbreviations and notations: we shorten G(no,io,jo) ^ and C{n(),io,jo) ^'^ ^- Recall 
the notation p(z,j) = no - K -«o| - |j -Jo I- So Fc(G) = {(z, j) : h{i, j) < p{i, j)}. 

Proposition 3. Let {i,j) E be a closed face of G and let {i',j') E J-' be such that i' is 
between i and iq and j' is between j and jq. Then {i',j') is also a closed face ofG. 

Proof. Without loss of generality, we may assume that io < i' < i and jo < j' < j. We have 
j') = p{hj) + — "^0 + (j — j')- Oil the other hand, we have h{i',j') < h{i,j) + {i — 
i') + (j - j')- As h{i,j) < p{i,j), we have h{i'J') < p{i',i'). □ 

So the faces of G form four "staircases" , as in figure 0] 
As a corollary, we deduce 

Proposition 4. G is connected. 

Proof. Clearly, if v and w lie on the same closed face of G, then v and w lie in the same 
connected component of G. But, by the previous proposition, we can travel from any closed 
face of G to (^o, jo) along a sequence of lattice-adjacent faces and lattice adjacent faces have 
a vertex in common. Thus, since every vertex of G lies on a closed face, every vertex must 
lie in the same connected component as the vertices of (io, jo)- So G is connected. □ 

It is also clear from the picture above that G is bounded by a single closed loop. Call 
this loop 5*. Divide 5*, as in figure |31 into four arcs by horizontal and vertical lines through 
(^05 Jo)- (In the diagram, 5* is in bold and the lines through (^o, jo) are dashed.) 
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Figure 4: Positions of Faces of a G(^n,i,j) 



Proposition 5. These lines divide S into four paths. Each of these paths contains an odd 
number of vertices. On each of these paths, the vertices alternate between vertices all of 
whose neighbors in Q are inside or on S (and thus in G) and vertices all of whose neighbors 
in Q are outside or on S. The ends of each path are of the latter type. 

Proof. We consider the section of S on which i > io and j > Jq, the other four sections 
are similar. For the faces {i,j) of Q within S (that is the closed faces of G) we have 
h{i,j) < riQ — {i — io) ^ (j Jo) where as for the faces which border S on the outside (that 
is, the open faces of G) we have h{i,j) = uq — {i — io) — (j — jo)- 

Our proof is by induction on the number of faces {i,j) of G for i > io and j > jo- In the 
base case, where there is only one such face, the path in question has one vertex and the 
result is trivial. When there is more then one such face, we can always find a face {i,j) such 
that (i + 1, j) and {i,j + 1) are outside S. For simplicity, we describe the situation where 
i > io and j > jo and leave the boundary cases to the reader. 

Set h = no — {i — io) — (j — jo)- We must have h{i + 1, j) = h{i,j + 1) = /i — 1 as they 
are open faces of G. We must have h{i,j) < h; as \h{i,j) — h{i + 1, j)| = 1 we must have 
h{i,j) = h — 2. If we then change h{i,j) to h, we remove (i, j) from within S. 

There are four cases, based on whether h{i — 1, j + 1) and h{i + 1, j — 1) = h or h — 2; 
figure El shows what happens in each case. The edges of S are drawn in bold, the faces of G 
are surrounded by thin boxes and the bipartite coloring of G is shown in black and white. 

□ 

We will use this proposition in the proof of proposition ^1 This proposition is also 
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useful in testing whether it is hkely that a graph can be reahzed through the method of 
crosses and wrenches. For example, the hexagonal graphs of section 6 of |Kuoj can not be 
so realized, because on their boundary there are six different boundary edges connecting 
vertices of degree two and a crosses and wrenches graph can only have four such edges. 
Examples 

3.6 Aztec Diamonds 

We begin with an example we discussed in the introduction: the Aztec Diamond graphs. 
In this case, h{i,j) = or — 1 when i + j mod 2 is or 1 respectively. In this case, every 
square of four values is 

0-1 -10 
or 

-10 -1 

so in every case we get a cross. 




The infinite graph Q is just the regular square grid. The edge labeling is exactly as described 
in section 11.21 The graphs G(^n,i,j) are the standard Aztec diamond graphs. 

3.7 Fortresses and Douglass' and Blum's Graphs 

In this section we show that Crosses and Wrenches graphs are capable of reproducing several 
previously studied families of graphs; the number of whose matchings are perfect powers or 
near perfect powers. Inside each face we put the number h{i,j). Most of these families 
of graphs appear to have first been described in print in |Ciucuj . which we will use as our 
reference. 

Fortress graphs are discussed in |Ciucuj , section four. To obtain a fortress graph, we take 

h{i,j) =0 i+j =0 mod 2, 

=1 i =0 mod 2, j = 1 mod 2 

= — 1 i =1 mod 2, j =0 mod 2. 
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Figure 6: A Portion of Q for Fortresses 



Every square of values is of the form 

1 0-1 10 



or 



-1 



-1 0' 1 ' -1 1 

so in every case we get a wrench, in the repeating pattern 



Joining up the wrenches, we get the infinite graph in figure IHl 

Once again, if we set all of the variables equal to 1, f{n,i,j) will count the matchings of 
fortresses. We must have f{h{i,j),i,j) = 1 and / must obey the defining recurrence. It is 
easy to check that both conditions are obeyed by 



f{2n + l,t,j) 
/(2n + l,«,j) 



5" 



2-5^ 



i = n mod 2 
j = n mod 2, 



a result of |Ciucuj . 
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Figure 7: An Example of a Fortress 

For example, figure [7| shows a fortress with 25 matchings, it is of the form 6*3,4 j where 
the center square is and i = mod 2, j = 1 mod 2. This fortress has 9 closed and 16 
open faces. 

Similarly, to get another family of graphs first studied by Chris Douglas (see |Ciucuj . 
section six), we take 

h{i,i)= i + i= mod 2, 
= 1 i + j = 1 mod 4 
=—1 i + j =—1 mod 4. 

We now sometimes will have 

10 -10 
or 

1 0-1 



which will produce a cross and sometimes will have 



-10 10 



or 

10-1 



which will produce a 
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Overall, we have the repeating pattern 



which produces the infinite graph in figure |H1 
Again, a quick induction verifies that 



/(2n,ii, j) 
f{2n + l,t,j) 
/(2n+l,^,j) 



i + j = 2n + 1 mod 4 
i + j = 2n — 1 mod 4. 



Similarly, to get the graphs considered by Matthew Blum (see |Ciucuj . section 8), take 

h{i,j) =0 i + j = mod 2 

= 1 i + j =1 mod 2, j = 0, 1 mod 4 
=—1 i+j=l mod 2, j =2, 3 mod 4. 

(Ciucu's methods require him to embed Q into the square grid, and he thus takes Q to be 
made of octagons and hexagons. For our purposes, it is more natural to apply Lemma [7| 
of Section 14.21 to obtain a grid of quadrilaterals and hexagons with the same number of 
matchings.) We get 



This gives us the infinite graph of figure IHl A quick induction then shows that 



/(3n,2, j) 
/(3n± 

f{3n±l,i,j) 



2-3 



n^±n 



3n± 1 


+ 


j 


3 


_2_ 


3n± 1 
3 


+ 


J 

_2_ 



1 mod 2 
mod 2. 
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Figure 9: A Portion of the Q for Blum's Graphs 
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3.8 Somos 4 and 5 



As described in the introduction, a combinatorial formula for the octahedron recurrence will 
give rise to a combinatorial interpretation of the Somos sequences. Specifically, if we set 
X = {{n, i,j) G C : —4 < 2n + i < 0} then the number of terms of /(O, 2k, 0) will obey the 
Somos-4 recurrence. Similarly, if X = {{n,i,j) G C : —5 < (Sn + i + 3j)/2 < 0} then the 
number of terms of /(O, 2k, 0) will obey the Somos-5 recurrence. 

Figures ITUl and ITT] show the infinite graph and the first several finite graphs corresponding 
to the Somos-4 recurrence. Figures IT^ and IT^ do the same for Somos-5. It may be hard to 
see the periodicity of the Q for (5, 1, 2); the tiles repeat as one travels five cells to the right, 
or one cell up and two to the right. 

4 Proof I: Urban Renewal 

In this section, we will prove the correctness of the crosses-and-wrenches algorithm. Our 
basic strategy is to hold the point (no, io,jo) where we are evaluating / fixed, while varying 
h. Our proof is by induction on the number of points in W fl Cno,io,jo- 

4.1 Infinite Completions 

In this section, we introduce an alternative way of viewing f{n,i,j) as counting the number 
of matchings of an infinite graph, subject to a "condition at infinity" to be described later. 
This means that a priori f{n,i,j) could contain an infinite number of terms, although in 
fact it will not, and describing the terms of f{n,i,j) requires an a priori infinite amount 
of information. On the other hand, this new method will no longer require the use of open 
faces and will remove the need to treat faces in the boundary as a special case in our proof. 
We will use this interpretation in our first proof of the Main Theorem. 

Let hhea. height function. Fix a triple (rio, "^o, jo) ^ ^, we will be interested in /(no, io, jo). 
Let G be the graph with open faces G(^no,io,jo)- Set 

h{i,j) = mm{h{i,j),p^rio,io,jo)ii^j))- 

It is easy to check that /i is a pseudo-height function. We can use the method of crosses 
and wrenches to associate an infinite graph Q to h. Note that if {n,i,j) is a face of G then 
h{i,i) = h{i,i). Thus, G is a subgraph of Q. Upon removing G, what remains looks like 
figured Let Gout = G\G. 

There is a unique matching (indicated in thick lines) of Gout given by taking the middle 
edge of every wrench. Call this matching Mout- Clearly, matchings of Q that coincide with 
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Figure 10: The Q for the Somos-4 Sequence 
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2 Matchings 

3 Matchings 



7 Matchings 



23 Matchings 

Figure 11: The First Four Nontrivial Somos-4 Graphs 

Mout on Gout are in bijection with matchings of G. Moreover, one can easily check that the 
face exponents associated to a given matching of G are the same as for the corresponding 
matching of including that the exponent of any face of Q that does not correspond to a 
face of G is 0. 

Thus, we could define f{no,io,jo) as the sum over matchings M of ^ such that M 
coincides with Mout on Gout of m(M). 

We claim that the following, less obvious result also holds: 

Proposition 6. Let M be a matching ofQ in which all hut finitely many vertices are matched 
to the same vertex as in Mout- (Mout is a matching of Gout, so for all but finitely many vertices 
of Q this makes sense.) Then M coincides with Mout everywhere on Gout- 
Proof. Suppose the opposite. All the vertices of G \ ^ come from wrenches. Each wrench 
has one of its vertices closer to (io, Jo) than the other; call these the near and far vertex of 
the wrench. Now, suppose the theorem is false. Of the finite number of vertices of ^ \ G not 
matched as in M, let v be the furthest from (io,io)- We derive a contradiction in the two 
possible cases. 

Case 1: f is a near vertex. Then the far vertex w in the same wrench as v is matched as 
in M . But w is matched with v in M, so w is matched with v and v is matched as in Mout- 

Case 2: w is a far vertex. But all of w's neighbors except the one it is matched to in Mout 
are farther from (io, Jo) than v is, so they are matched as in Mout and can not be matched 
to V. So is matched as in Mout- n 

As result, we can give another statement of the main theorem. 
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Figure 12: The Q for the Somos-5 Sequence 
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5 Matchings 



1 1 Matchings 



Figure 13: The First Four Nontrivial Somos-5 Graphs 
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Figure 15: Splitting a Vertex 

The Main Theorem (Using Infinite Completions). Call a matching M of Q an ac- 
ceptable matching if M coincides with Mout for all hut finitely many edges. Then 

/K,^o,jo)= X] m(M). 

M acceptable 

4.2 Some Easy Lemmas 

Proposition 7. Consider a graph with open faces G. Replace a vertex G with two unweighted 
edges as in figure (731 to form a new G' . 
Then M{G) = M{G'). 

Proof. Any matching of the old graph can be uniquely transformed to a matching of the new 
graph by adding one of the two new edges. The additional edge in the matching contributes 
a factor of 1 in the product of the edge variables. The two faces adjacent to the new edges 
have one more used edge and one more unused edge than the corresponding faces of the 
previous matching, so they have the same exponent. The other faces are unaffected. □ 

The Urban Renewal Theorem. Suppose a graph with open faces G contains the sub-graph 
with open faces shown in the left of fiaureURwith the indicated edge weights and face weights. 
(The face weights are in uppercase, the edge weights in lower case.) Create a new graph G' 
by replacing this with the graph with open faces in the right of fiaure\l(A where 

X' = {adWY + hcVZ)/X. 

Then m{G') = m{G). 

This theorem is related to Lemma 2.5 of |Ciucuj . However, because he doesn't use face 
weights, Ciucu's theorem involves slightly different replacements and yields a relation of the 
form m{G) = (ac + bd)m{G'). 

Proof. Write m{G) = niQ + nii + m2 where mo consists of the terms corresponding to 
matchings where none of the edges a, b, c or d is used, mi consists of the terms using one 
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Figure 16: Applying Urban Renewal 
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Figure 17: The Bijection in the mi Case 



^ D 

Figure 18: The One-to-Two Map in the mo Case 

such edge and m2 consists of the terms using two. Write m{G') = + m'^ + m'2 similarly. 
We will show that mi = m'^, mo = and m2 = itlq. 

We first give a bijection between the matchings counted by mi and those counted by m'^ 
as shown in figure El 

In each case, it is easy to check that the two matchings paired off have exponent on X 
and X' and raise all other variables to the same exponent. 

Next, as shown in figure ITHl we give a bijection that associates to each matching counted 
by mo two of the type counted by m'2. 

The first transformation in figure ITHl increases the exponents of W and Y by 1, adds the 
edges a and d and changes the contribution of the center face from X to X'^^ . Similarly, 
the second transformation increases the exponents of V and Z by 1, adds the edges h and c 
and changes the contribution of the center face from X to X'~^ . So we get from mo to 
by multiplying by 

adWYX''^ + hcVWX'^^ 

X ■ 

But, as X' = {adWY + bcVZ)/X, this is 1. 

Finally, we show 1712 = mp. This is similar to the preceding paragraph. This time, we 
pair off every two matchings in 1712 with one matching in mg, as shown in figure IT^ 
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□ 

□ 



Figure 19: The Two-to-One Map in the m2 Case 



In the first pairing of figure UHl we delete edges a and d, reduce the exponents of W and 
F by 1 and replace by X'. In the second, we delete edges b and c, reduce the exponents 
of V and Z by 1 and replace X^^ by X'. Thus, we are multiplying by 

X' 



adWYX-^ + hcVWX-^ 



which is 1 as before. 



□ 



4.3 Proof of the Main Theorem 

Let (no, jo) £ W and let h and h be as in subsection 14. 1| let Q be the graph with open faces 
associated to h. Our proof is by induction on j{p{i, j) — h{i,j)) = 2#(W fl C(no,io,jo))- 

li'^i j{p{i, j) — h{i, j)) = then h{i, j) =p{i,j) soQ is the graph shown in figure QUI which 
has only the indicated admissible matching. The matching polynomial of this matching is 
x{io,jo)- We also have in this case that /i(«0) Jo) = p(^O)Jo) = ''^o so we have f{no,io,jo) = 

If J2ijPi'^^j) ~ f^ihj) > then, among the (z,j) for which h{i,j) < p{i,j) there is (at 
least) one with h{i,j) minimal, let this be We claim that h{i ± 1, j) = h{i,j ± 1) = 

h{i,j) + 1. We deal with the case of h{i + 1, j), the other three cases are similar. We have 
h{i + 1, j) = h{i,j) ± 1 as is a pseudo-height function. If h{i + 1, j) = h{i,j) — 1, then 

+ 1> j) < Hhj) < Pihj) - l<pii + IJ), 

contradicting the minimality of 

Thus, we have shown h{i±l,j) = h{i,j±l) = h{i,j) + l. So the face associated to h{i,j) 
is a square. Apply the Urban Renewal Theorem to this square to create a new graph with 
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Figure 20: The Base Case, with its Unique Admissible Matching 



the same matching polynomiaL Then apply Lemma [7| to this graph wherever possible. The 
graph thus created will be the same as replacing any crosses whose vertex lies on the face 
by a wrench with one vertex on and vice versa. Also, the edge weights adjacent 
to are interchanged with their diametric opposites and the weight x{i,j) is replaced by 
a certain binomial. 

The graph thus produced is precisely the graph produced by the function h' where 
h'{i,j) = h{i,j) + 2 and h' = h everywhere else. By induction, the matching polynomial of 
this graph is f{nQ,iQ, jo) for the initial conditions h'. The / for the initial conditions h is 
given by replacing x{i,j) by the same binomial as before. Thus, /(rio, io,jo) is precisely given 
by the matching polynomial of the graph for h' with x{i,j) replaced by the said binomial, 
which by the Urban Renewal Theorem is precisely the matching polynomial of the graph for 
h. □ 

5 Proof II: Condensation 

In this section, we give a second proof of the main theorem, similar to the proof of |Kuoj . In 
order to avoid lengthy tedious verifications, we will perform this proof with the face variables 
set equal to 1 throughout this section. 
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5.1 The Condensation Theorem 

We first give some graph theoretic notations. 

Let G be a graph and S C V{G). Let d{S) denote the elements of S that border vertices 
of G not in 5*. If G is bipartite, colored black and white, let 6{S) denote the number of black 
vertices of 5* minus the number of white vertices. While S{S) is only defined up to sign, we 
adopt the implicit assumption that, if Si and 5*2 are both subsets of V{G), then S{Si) and 
6(82) are computed from the same coloring of G. Let 5^(5") be the subgraph of G induced 
by S. We abuse notation by writing m{S) to mean m((7(S')). (Recall that all face variables 
have been set equal to 1, so m{S) is a polynomial in the edge variables of g{S).) 

The following theorem is essentially due to Kuo and proven in jKuo] . (Kuo's paper proves 
this relation for many specific families of graphs but never states the general theorem. The 
below is my attempt to assimilate all of Kuo's cases into a general framework which will also 
encompass the results of this paper.) 

Kuo's Condensation Theorem. Let G be a bipartite planar graph and let the vertices of 
G be partitioned into nine sets 

V{G) = CUNUNEUEUSEUSUSWUWU NW. 

Assume that only vertices in the sets joined in figure Wl\ border each other. (Vertices may 
also border other vertices in the same set.) Assume further that 

5(NE) = 5(SW) = 1, 5(SE) = 5(NW) = -1 

and 

5{C) = 5(N) = 5(E) = 5{S) = 5(W) = 

Finally, assume that d(NE) and (9(SW) are entirely black and 5(NW) and (9(SE) are 
entirely white. 
Then we have 

m{G)m{C) = m(N U NE U NW U C)m(S U SE U SW U C)m(E)m(W) + 

m(E U NE U SE U C)m(W U NW U SW U C)m(N)m(S). 

We write X to denote the collection of sets {N, NE, E, SE, S, SW, W, NW}. 

Proof. Let G always denote a graph obeying the above hypotheses. We define a northern 
join to be a set 7 of edges of G of any of the following types 

1. A path with one endpoint in NE, the other endpoint in NW and the intermediate 
vertices in C. 
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2. A single edge with one endpoint in NE and the other in NW. 

3. A pair of edges, one joining NE to N and the other joining NW to N. 
We define an eastern, southern or western join analogously. 

Proposition 8. Let M be a multi-set of edges of G such that every member of C lies on 
two edges of M and every other vertex lies on a single edge of M. Then we can write M 
uniquely as a vertex disjoint union 

M = 7i U 72 U Mc U y Mq 

where 

1. Either 71 is a northern join and 72 is a southern join, or 71 is a eastern join and 72 
is a western join. 

2. Mc is a disjoint union of cycles, entirely contained in C. (We count a doubled edge 
as a cycle of length 2.) 

3. Mq is a vertex disjoint set of edges contained entirely within q. 

Proof. Since every vertex is adjacent to either one or two elements of M, M can be written 
uniquely as a disjoint union of cycles and paths. Moreover, the vertices of C are on cycles 
or are inner vertices of paths and the vertices of V{G) \ C are the ends of paths. 

Let Mc be the cycles of M, Mc is vertex disjoint from the other edges of M. After 
deleting Mc, the remaining graph is a disjoint union of paths, each of whose endpoints lie 
in one of the q e X and all of whose interior vertices must he in C. Let Mq be the isolated 
edges connecting one point of q to another. After deleting all of the Mq, what remains 
are paths as before which either have at least one interior point or whose endpoints lie in 
different q G X. Call the set of remaining edges A. 

Let q G X. We claim A fl q C d{q). This is because, if f G q \ d{q), then there is a 
single edge vw of M containing v. We have w G q by the definition of 9(q). Then vw is 
also the only edge containing w and thus vw is an isolated edge of M lying entirely in q. So 
V G Mq and q ^ A. 

We know Mne contains equally many black as white vertices as it is a disjoint union of 
edges. Thus A fl NE has one more black vertex than white, as we assumed 5(NE) = L 
However, we just showed A fl NE C 5(NE), and 9(NE) is entirely black. So there is a 
unique vertex i>ne in AflNE. Similarly, there are unique vertices vjsrw, vse and vsw, with 
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vne and t'sw black and ^nw and vse white. Also by the same logic, #(A fl q) is even for 
q G {C, NE, NW, SE, SW}. 

Now, the Vq must be endpoints of paths. The a priori possible paths are the following, 
and their rotations and reflections: 

1. A path from vjsiE to t'NW through C. (There are four possibilities in this symmetry 
class. 

2. A path from to vsw through C. (There are two possibilities in this symmetry 
class.) 

3. A single edge from v-^e to vnw- (There are four possibilities in this symmetry class.) 

4. A single edge from to A^. (There are eight possibilities in this symmetry class.) 

As #(A n N) is even, there must either be or 2 paths of type (4) ending in N, and 
similarly for E, S and W. This means we can not have a path joining wne to Vsw- If such 
a path existed, v-^w could not be joined to or W, as that would create one path ending 
in that set. Also, t'NW could not be joined to vse, as the paths i'neI'sw and i'nwI'se would 
cross. Similarly, we may not have a path joining vjsiw to vse- So the paths of type 2 do not 
occur. 

It is now easy to see that A must decompose either as the union of a northern and a 
southern join or as the union of an eastern and a western join. □ 

Note that, if 7^ (i = 1 or 2) passes through C, it will contain an odd number of edges, as 
its endpoints are of opposite colors. 

We now prove the result. Each side of the equation is a sum of products of edge variables, 
and each product corresponds to an M meeting the conditions of the lemma; we will abuse 
notation and call this product M as well. We will show each M occurs with the same 
coefficient on both sides. 

Specifically, let k denote the number of cycles in Mc of length greater than 2. We claim 
that M appears on each side of the equation with coefficient 2*^. Even more specifically, we 
claim that m(G)m(C) contains M with coefficient 2'^. We also claim that, if 71 is a northern 
join and 72 a southern join, then m(N U NE U NW U C)m(S U SE U SW U C)m(E)m(W) 
contains M with a coefficient of 2^= and m(EUNEUSEUC)m(WUNWUSWUC)m(N)m(S) 
does not contain M. If 71 is an eastern and 72 a western join then we claim the reverse holds. 

Consider first the task of determining how many times M appears in m{G)m{C). This 
is the same as the number of ways to decompose M as M = M{G) U M(C) with M{G) and 
M(C) matchings of G and C respectively. 
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Clearly, the edges with endpoints outside C must lie in M{G). This means all the edges 
of Mq, q G X, the edges of 7^ if 7^ does not pass through C and the final edges of 7^ if 7^ does 
pass through C. This forces the allocation of all the edges of 7^ because, when two edges of 
M share a vertex, one edge must go in M{G) and one in M(C). This forcing is consistent 
because 7^ has an odd number of edges and the end edges must both lie in M{G). Finally, 
we must allocate the edges of Mq- All the cycles of Mq are of even length as G is bipartite. 
In a cycle of length greater than 2, we may arbitrarily choose which half of its edges came 
from M{G) and which from M(C). Thus, we make 2^ choices. 

Next, consider the coefficient of M in m(N U NE U NW U C)m(S U SE U SW U 
C)m(E)m(W). We must similarly write 

M = M(N U NE U NW U C) U M(S U SE U SW U C) U M(E) U M(W). 

We claim that if 71 is an eastern join and 72 a western, this coefficient is 0. There are three 
cases. 

Ccise 1. 7i is a path joining NE to SE through C. 

Then every edge of 71 must lie in M(N U NE U NW U C) or in M(S U SE U SW U C) 
and which of these two it lies in must alternate. There are an odd number of edges in 71, so 
the final edges must lie in the same one of these two sets. But any edge with an endpoint in 
NE must lie in M(N U NE U NW U C) and any edge with an endpoint in SE must lie in 
M(S U SE U SW U C) U M(E). 

Case 2. 71 is a single edge with endpoints in NE and SE. 

But none of the M(q)'s can contain such an edge. =^<^ 

Case 3. 71 is a union of two disjoint edges, one connecting NE to E and one connecting 
SE to E. 

But neither of these edge types can lie in any of the four M(q)'s. ^<^= 
Finally, we show that if 71 is a northern join and 72 a southern join that M has coefficient 
2^ in m(N U NE U NW U C)m(S U SE U SW U C)m(E)m(W). As before, the edges of Mq, 
q,e X, are uniquely allocated to one of M(N U NE U NW U C), M(S U SE U SW U C), 
M(E) and M(W). If 7^ does not pass through C, its edges are also immediately forced. If 
7j does pass through C, its final edges are forced, thus forcing the others and again we have 
consistency as the two final edges are in the same set and 7^ has an odd number of edges. 
Finally, each cycle of Mq of length greater than 2 can be allocated in two ways. 

The case of m(EUNEUSEUC)m(WUNWuSWuC)m(N)m(S) is exactly analogous. 

□ 



41 



5.2 The Decomposition of G 



Let G be a graph arising from the method of crosses and wrenches for some height function 
h and some (nc^cJo)- The purpose of this section is to describe a decomposition of the 
vertices of G into nine disjoint sets. In the next section, we will show these sets obey the 
hypotheses of the previous section. For simplicity, we assume that (no — 2, zq, jo) ^ ^■ 
Note that we have 



y{G{no-2,io,jo)) - ^(^(no-l,io+l,io))) l~l ^(CCno-l.io-lJ))) 
= y{G{no-l,io,jQ-l)) l~l ^(<^(no-l,iojo-l)) 

and 

^(G(no,jojo)) = ^(^(no-l,io+l,io)) U ^^(^(no-i^io-lj)) 

(Recall that V{G) denotes the vertices of G.) 

Thus, the four graphs G(no — 1, ± 1, jo) and Gin^ — 1, zq, jo ± 1) intersect as shown in 
the Venn diagram in Figure |22l where the unlabelled cells indicate empty intersections. We 
decompose \/(G(no, ^o, jo)) into the nine sets C, N, NE, E, SE, S, SW, W and NW as 
shown in the figure. 

We now verify that the sets of vertices C, N, E, W, S, NE, NW, SE and SW obey 
the hypotheses of the Condensation Theorem. We also show that m(E) = a(io + — 1, jo), 
m(W) = c(2o-?^o + l, jo), m(N) = 6(io, jo + ^^0 - 1), "^(S) = (i(io, jo-^o + !)• These clearly 
establish the Main Theorem. 

The following lemma is clear. 

Proposition 9. Abbreviate p (no, io, jo) {h j) byp{i,j). Then 

{hi) e i^c(G'(„(,-2,»ojo)) ^// Khj) < Pihj) - 2 

) iff h{i,j) = p{i,j) -2 and i > io 

(no-2,io,jo) )) iff Khi) = P{hj) - 2 and j > jo 

(i, j) e F^{G 

(no-l,io-l,jo)) \ -^c(C(no-2,io,jo) ) ^// h{i,j)=p{i,j)-2and i < io 

G Fc{G 

(no-l,io,io-l) ) \ Fc{G 

(no-2,io,jo) ) iff h{iJ) = P{iJ) - 2 and j < jo 

(We are using our simplifying assumption that (^o, jo) ^ Fc{G(no-2,io,jo)) avoid a messy 
statement for this face.) 
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Figure 23: The Geometry and Coloring of E 



Proposition 10. We have m(E) = a(io + no — I, jo), m(W) = c{io — no + 1, jo); '^(N) = 
K'^o,3o + no-l) andmiS) = d{io,jo-no + l) and6{E) = 6(W) = 6{N) = 6{S) = 0. Also, all 
of the vertices of E that border vertices of NE are the same color. Calling this color white, 
the vertices of N bordering NE, those of W bordering SW and those of S bordering SW 
are white. The vertices of E bordering SE, those of S bordering SE, those of N bordering 
SW and those of W bordering SW are black. 

Proof. We do the case of E, the others are similar. Let v & £ and let j) be a closed 
face of bordering v (by the definition of E, such a face exists); by Lemma 01 

h{i,j) = p{i,j) — 2. Now, is not a closed face of G(^no-i,ioM+i)- Thus, by Lemma El 
j < Jo- But, similarly, j > jo so j = jo- Moreover, the faces that v borders which are not of 
the form {i,jo) must not lie in G at all, so we have h{i,j ± 1) = h{i,j) + 1. 

Putting this data into the Crosses and Wrenches algorithm, we have that E looks like 
figure 1221 where the number of hexagons and the length of the dangling path can vary. The 
highlighted vertices are those that border NE and SE. 

Clearly, E has a unique matching M, shown in bold. Let i be the largest value for 
which {i,jo) is a closed face of G, the only weighted edge in M is the edge separating 
(i, jo) and (i, jo + 1). As h{i,jo) = {no - 1) + {io + I - i) - (jo - jo) = no + io - i and 
h{i + 1, jo) = h{i,jo) + 1 = no — io + * + 1 (as otherwise (i + 1, jo) we be a closed face) we 
see that this is the edge 0(^0 + '^o ~ 1; jo)- As E has a matching, 5(E) = 0. 

All that is left to show is that the vertices of N that border NE and the vertices of E 
that do are the same color. It is equivalent to show that the vertex of b{i,j + n — 1) and 
that of d{i + n — 1, j) closer to NE are the same color. But the path joining them is the one 
we proved in Proposition El had an odd number of edges. □ 

Proposition 11. Let v G NE. Then exactly one of the following is true. 



1. V lies in a wrench oriented NE-SW (i.e., this orientation: r ), the other end of which 
is also in NE. 

or 

2. V lies on the outer perimeter of G^no^io.jo)- 
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Proof. Let {i + 1, j), + 1), {i + 1, j + 1) be the faces surrounding v {y may only be 

adjacent to three of them) and let I be the set of these four faces.. By Lemma IHl either i or 
i + 1 > iQ so i > iQ. Similarly, j > jo- There are several cases. 

First, suppose h{i,j) = p{i,j) — 2. Then h{i',j') > p{i',j') — 2 for G /. If 

h{i',j') = p{i',j') — 2 for all G /, then v lies on a NE-SW wrench both ends of which 

are in NE. If some h{i',j') > p{i',j') — 2 then is not a closed face of G(^no,io,jo) 

one can check that in every case v is adjacent to so v is on the outer perimeter. 

Now, suppose that h{i,j) < p{i,j) — 2. Then G Fc{G(no-2,io,jo))- The only way 

that V could be in E then is if v were at the NE end of a NE-SW wrench so that it didn't 
border Then {i + does border v, so h{i + 1, j) > p{i + 1, j) — 2. The only 

way this is possible is if h{i,j) = p{i,j) — 4, h{i + = h{i,j + 1) = p{i,j) — 3 and 
h{i + 1, j + 1) = p{i,j) - 2 = p{i + l,j + 1). Then {i + 1, j + 1) ^ FciG(no,io,jo)) and v lies 
on {i + 1, j + 1), so again v is on the perimeter. 

Finally, suppose that h{i,j) > p{i,j) — 2. Then h{i',j') > p{i',j') — 2 for all G /. 

But then none of the are closed faces of G(^no,io,jo) v is not a vertex of G(^no,io,jo)^ ^ 

contradiction. □ 

Proposition 12. NE and SW have one more black vertex than white, and vice versa for 
SE and NW. 

Proof. We discuss the case of NE, the others are similar. 

Pair off all of the vertices of E for which the first case of the previous proposition holds 
with the other end of their wrench. This removes equal numbers of black and white vertices. 
What remains is a path, which we showed in Proposition 1101 has each end black. □ 

Proposition 13. (9(NE) and (9(SW) are entirely black and vice versa for (9(SE) and 
(9(NW). Moreover, only the adjacencies permitted by the hypotheses of the Condensation 
Theorem can occur. 

Proof. Besides the endpoints of the path in the previous proposition, 9(NE) must be made 
up entirely of vertices from the wrenches in the previous proposition. However, in fact, only 
the south-west vertex of a wrench will be able to lie in (9(NE), and these are all the same 
color. 

This deals with the coloration of the boundaries. This immediately means that NE 
can not border SW. We saw in Lemma ^1 that NE only borders N and E. These and 
the symmetrically equivalent restrictions gives the required restrictions on what may border 
what. □ 

Proposition 14. C has equally many black and white vertices. 
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Proof. C = G(„o_2,io,jo). By induction on uq, we can assume that we already know that 
/(no — 2,io, Jo) — '>Ti{C). In particular, C has at least one perfect matching. □ 

We have now checked all of the conditions to apply the Kuo Condensation Theorem. 



5.3 A Consequence of the Condensation Proof 

The Condensation proof of the main theorem, although it involves more special cases, seems 
to have at least one consequence which does not follow from the Urban Renewal proof. 

Proposition 15. In the relation 

f(n, i, j)f(n - 2, i, j) = a(i + n - 1, j)c{i - n + 1, j)f{n + l)f(n - 1, i, j - 1) + 

Kh j + n- l)d{i, j -n + l)f{n - 1, i + 1, j)f{n - 1, i - 1, j) 

each monomial appearing on the right hand side appears in exactly one of the two terms on 
the right, and appears with the coefficient 2^^ for some k. 

Proof. This just says that every combination of edge M as in the proof of the Kuo Conden- 
sation Theorem occurs either only in 

m(N U NE U NW U C)m(S U SE U SW U C)m(E)m{W) 

or only in 

m(E U NE U SE U C)m(W U NW U SW U C)m(N)m(S) 

and appears with coefficient 2'^ for some k; this follow from the proof of the Kuo Condensation 
Theorem. □ 



6 Some Final Observations 

6.1 Setting Face or Edge Variables to One 

Proposition 16. If we replace all of the face variables in f{n,i,j) by 1, every coefficient in 
f{n,i,j) will still be 1. 

Proof. We are being asked to show that we can recover a matching of Gn,i,i (henceforth, G) 
from knowing which weighted edges were use in the matching. Since the unweighted edges 
of G are vertex disjoint, knowing which weighted edges are used determines all the edges 
used, and hence the matching. □ 
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The corresponding version for faces is true but more difficult, and its proof has more 
interesting consequences. 

Proposition 17. If we replace all of the edge variables in f{n,i,j) by 1, every coefficient in 
f{n, i,j) will still be 1. 

Proof. We are now being asked to show that we can recover a matching M of G from just 
knowing the face exponents. By the previous result, it is enough to show that we can 
determine whether or not a weighted edge appears in M. Recall the notation of Section 12.21 
for / a face of e(/) is the exponent oi Xf and, for e an edge of 5{e) is the exponent of 
e. 

We give a formula for 5e in terms of certain e{f)'s in the case where e is of the form 
{a,i,j), similar formulas for 6{b,i,j), 6{c,i,j) and S{d,i,j) can be found and will be stated 
later in the section. My thanks to Gabriel Carrol for suggesting this simple formulation and 
proof of the lemma (in a somewhat different context.) 

Proposition 18. Assume that Gn,i,j contains the edge (a, iojjo)- IiT' o-ny matching of G, we 
have 

(^(a,«o, jo) = - ^ e{n,i,j). 

{n,i,j)£l 
n+i+j<io+jo+l 
n+j-j<io-io+l 

Proof. Let t denote a formal variable independent of all other variables. 
Let / be defined by 

f{n,i,j)f{n - 2,i,j) = a{i + n - l,j)c{i-n + lj)f{n - l,i,j + l)f{n - l,i,j - 1) + 

b{i,j + n- l)d{i,j -n + l)f{n - l,i + l,j)f{n- l,i- 1, j) 

for the initial conditions x{i,j) on I. 
Define / by 

/(n, i, j) = tf{n, i,j) n + i+j<io+jo + l and n + i-j<io-jo + l 
fin,ij) = f{n,i,j) otherwise. 

Define a(io) Jo) = tcL{hj) ^ind a{i,j) = a{i,j) otherwise. Then we claim that 

fin,i,j)f{n - 2,i,j) = a{i + n - l,j)c{i-n + 1, j)/(n - + l)f{n - 1,2, j - 1) + 

b{ij + n- l)d{i,j -n + l)f{n -l,i + 1, j)/(n -l,i- 1, j). 

This may be checked by checking all nine cases for how n + i + j and n + i — j compare 
to io + jo + 1 and zq — jo + 1 , and noting that n + i + j = io + jo + 1 and n + i—j = iq— Jq + I 
imply that i + n — 1 = i^ and j = jo. 
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Now, let us consider the coefficient of t in any term of f{n, i, j) . Since (a, io, jo) appears in 
Gn,i,j, we must have n-\-i+j > + 1 and n+i—j > io—jo + l. Thus, f{n, = f{n, 

and no t appears. On the other hand, / can also be obtained from / by substituting tx{i,j) 
for x{i, j) for every (n, i, j) G X such that n + i + j < io + Jo + 1 and n + i—j < io — jo + 1 
and substituting ta{io,jo) for a(io, Jo)- So we see that the exponent of t is 

S{a,hj)+ ^ ^{n,i,j). 

{n,i,j)eJ 
n+i+j<io+jo+l 
n+i-j<io-jo+i 

So this sum is zero and we are done. □ 

So the face exponents determine the edge exponents and (by the previous Theorem) 
determine the matching. □ 

One can find three more formulas for S{a, i,j) and four corresponding formulas for each 
of 5{b,i,j), 5{c,i,j) and S{d,i,j). Specifically, let 

P = {{n,i,j) eT -.n + i+j <io+jo + l, n + i - j < io - jo + 1} 

Q = {in,i,j) el -.n + i+j <io+jo + l, n + i-j >io-jo + l} 

R = {{n,i,j)eT:n + i+j>io+jo + l,n + i-j>io-jo + l} 

S = {{n,ij) el -.n + i+j >io+jo + l: n + i- j <io- jo + l} 



Then 

^(a,i,j) = - ^ e{n,i,j)= ^ e(n,i,j) = l- e{n,i,j)= ^ e{n,i,j) 

{n,i,j)&P {n,i,j)eQ {n,i,j)eR {n,i,j)eS 

The same result holds for 

1. 5{b, i, j), if we define P, Q, R and S by comparing n ± i + j + 1 to ±io + jo- 

2. 5{c, i, j), if we define P, Q, R and S by comparing n — i ± j + 1 to —io ± jo- 

3. d{d, i,j), if we define P, Q, R and S by comparing n ± i — j + 1 to ±io — jo- 

A geometrical description can be given for these rules. Let e be a weighted edge of G. 
Draw four paths in G by starting at e and alternately turning left and right. (There are two 
ways to leave e and two choices for which turn to make first.) These divide the faces of G 
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into four regions P, Q, R and S. If the edge e is present, then the sums X](n i j)GP ^(^' "^5^)' 
T.{n,i,j)eQ4^^hj), E(n,ij)Gi?e(^,^,j) and E(„,i,j)g5 ^l"-, ^, j) will yield a 1, two -I's and a 0. 
If e is absent, one of these sums will be 1 and the other three will be 0. 

Interestingly, if one follows the operation of the previous where e is an unweighted edge, 
one gets precisely the same conclusions except that the interpretation of the results is changes: 
getting the sums (1,0,0,0) (in some order) now means that e is present and (0,1,-1,1) 
means e is absent. 

6.2 Height Functions for Crosses and Wrenches Graphs 

Let G be a bipartite connected planar graph with a fixed labeling of its vertices as black and 
white and let w be a function assigning a real number to each edge of G. We define a Propp 
height to be a function H from the faces of G to the real numbers with the following property: 
Let / and /' be faces of G separated by an edge e such that, when one stands on / and looks 
toward /', the white vertex of e is on the left. Then H{f) — H{f') = w{e) + (0 or — 1). 
Define two Propp heights H and H' to be equivalent if H{f) = H'{f) + c for some constant c 
independent of /. (The term Propp height is my own, to distinguish them from the heights 
which occur through out this paper.) 

In section 3 of |Proppl| , Propp essentially proved 

The Correspondence Between Heights and Matchings. For G as above, one can 
find w such that there is a bisection between equivalence classes of Propp heights for w and 
matchings of G. The bisection is such that, if H corresponds to M, the edge e occurs in M 
iff H{f) = H{f') + (w(e) — 1) where e separates f and f. 

These Propp heights have been very useful in studying statistical properties of random 
matchings, see e.g. |CEPj . For crosses and wrenches graphs, the Propp height has an 
extremely simple description. 

Proposition 19. Let Q be a crosses and wrenches graph. We may take w in the previous 
theorem to be given by w{e) = 1/4, if e is a horizontal or vertical edge and w{e) = 1/2, if e 
is a diagonal edge. 

Proof. We just must verify that, for every vertex v of Q, J2e5v'^(^) ~ where the sum is 
over e all incident on v. This is clear. □ 

In this context, the condition at infinity for infinite completions can be stated as H{i,j) = 
{\i\ — |j|)/4 for all but finitely many {i,j). 



49 



6.3 Generating Random Matchings 

Let G be a crosses and wrenches graph. Suppose that we want to randomly generate a 
matching of G, with our sample drawn from all matchings of G with uniform probability. 
Such sampling has produced intriguing results and conjectures in the cases of Aztec diamonds 
and fortresses. 

It is possible to do so in time 0{\G\'^), where |G| can be any of \F{G)\, \E{G)\ or \V{G)\ 
as all of these only differ by a constant factor for G a crosses and wrenches graph. More 
specifically, if G = G(^no,io,jo)y ^o so in time 0{\C[no,io,jo) nW|). (In most practical 

cases, this is closer to |G|^/^. 

We give only a quick sketch; the method is a simple adaptation of the methods of |Propp2| . 
Let G = G(^no,io,jo) arise from a height function h and let F = F{G). Throughout our 
algorithm, x will denote a function F ^ M. Our algorithm takes as input h and x and returns 
a random matching, where the probability of a matching M being returned is proportional 
to UfeF^ifY^^^- We describe our output as a list of the weighted edges used in M. At the 
beginning of the algorithm, we take x{f) = 1 for all /. 

Step 1: If /?.(zo) Jo) = "^O; return {} and halt. 

Step 2: Find G F such that h{i,j) = h{i ± 1, j ih 1) — 1 for all four choices of the 

±. 

Step 3: Set x'{i,j) = {x{i + l, j)x{i — l, j) + x{i, j + l)x{i, j — l))/x{i, j) and x'{f) = x{f) 
for all other / e F. Set h'{i,j) = h{i,j) + 2 and h'{f) = h{f) for all other /. 

Step 4: Run the algorithm with input h' and x', let the output be M' . 

Step 5: For shorthand, set a = {i + 1 + h{i,j),j,a), b = + 1 + h{i,j),b), c = 
{i — 1 — h{i,j),j, c) and d = — 1 — h{i,j), d). If M' contains one of {a, 6, c, d}, return M 
and halt. If M' contains two of {a, 6, c, ci}, return M' \ {a, 6, c, d} and halt. If M' contains 
none of {a, 6, c, rf}, return M U {a, c} with probability j + l)x{i,i — l)/x{i,j)x'{i,i) and 
return M U {6, c?} with probability x{i + — l,j)/x{i,j)x'{i,j). All other cases are 

impossible. 

In order to not do a time consuming search at step 2, one can keep a reverse look-up 
table which, given h, record all with h{i,j) = h and is updated in step 3 when h is 
replaced by h. Note that x{f) will always be integer valued. 
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